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Abstract 

Weak interaction rates and the associated energy losses for massive nuclei in the /p-shell are calcu- 
lated in stellar matter using a modified form of proton-neutron quasiparticle RPA model with separable 
Gamow- Teller forces. A total of 209 nuclei with mass numbers ranging from A = 61 to 80 are consid- 
ered here for the calculation of weak rates. These also include many neutron-rich nuclei which play a 
key role in the dynamics of the core collapse calculations. The stellar weak rates are calculated over a 
wide range of densities (10< pYe (gcm"^) < lO") and temperatures (lO'' < T(K) < 30 x 10^). This 
is our third paper in a series starting with the calculation of stellar weak rates for nuclei in the sd-shell 
[ATOMIC DATA AND NUCLEAR DATA TABLES 71, 149 (1999)]. The calculated rates take into 
consideration the latest experimental energy levels and ft value compilations. The effect of particle 
emission processes is taken into account and the energies and probabilities of these particle emission 
processes are also calculated in stellar environment. Our results are presented here on an abbreviated 
scale of temperature and density. 
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1 INTRODUCTION 



the shock energy 



Gravitational collapse of the core of massive stars 
(M > IOMq) is triggered by electron capture on 
free protons and on nuclei at high temperatures 
and densities. These stars then end up either as a 
neutron star or a black hole by supernova explosion 
(Type II). Electron capture during the late stages 
of stellar evolution was discussed in detail by [Q. 
The weak interaction certainly plays a key role in 
the dynamics of the core collapse. Hundreds, even 
thousands of nuclear reactions immediately before 
and after the collapse are relevant, and many of 
these concern unstable nuclei that are difficult to 
study experimentally. Therefore detailed, reliable 
theory is needed in order to understand these vio- 
lent processes. 

Because the final outcome of the explosion de- 
pends so sensitively on a variety of physical inputs 
at the beginning of each stage of the entire process 
(i.e., collapse, shock formation, and shock propa- 
gation), it is desirable to calculate the presuper- 
nova stellar structure with the best possible physi- 
cal data and input currently available. The energy 
budget would be balanced in favor of an explosion 
by a smaller precollapse iron core mass and a lower 
entropy. Since the hydrodynamic shock, believed to 
be responsible for the explosion, forms at the edge 
of the homologous core (the inner part of the core 
which remains in acoustic communication through- 
out core collapse and bounce), a smaller iron core 
size means that the energy lost by the shock in pho- 
todisintegrating iron nuclei in the overlying matter 
would be correspondingly smaller As less en- 
ergy is stored in nuclear excited states in the lower 
entropy environment of the presupernova as well as 
the collapsing core, the collapse can proceed to a 
higher density and generate a stronger bounce and 
form a more energetic shock wave [Q. In addition, 
at smaller entropy the abundance of free protons 
remains less. Since protons are the main sink for 
electrons through the electron capture reactions, 
a smaller initial entropy even at the presupernova 
stage means a lower overall electron capture rate 
and therefore a higher value of Yg (the number of 
electrons per baryon) at the time of bounce. This 
in turn makes the mass Mhc of the unshocked in- 
ner core (the homologous core) larger. As a result, 



c.r/;/^(i;/-i;,), 

is larger for larger final lepton fraction Y^f^ where 
Rhc is the homologous core radius and Yei is the 
initial electron fraction. 

The nuclear weak-interaction mediated reac- 
tions, e.g., P decay and electron capture , not only 
lead to a change in the neutron-to-proton ratio in 
the stellar core material but because of the removal 
of the energy by neutrinos produced in the reac- 
tions, they cool the core to a lower entropy state. It 
is therefore important to follow the evolution of the 
stellar core during its late stages of hydrostatic nu- 
clear burning with a sufficiently detailed nuclear re- 
action network that includes these weak-interaction 
mediated reactions. 

In the aftermath of the collapse, r-process nu- 
cleosynthesis (theory of which was put forward by 
^) hinges on the waiting-point nuclides, where 
the P decay rate is long enough for substantial 
amounts of material to build up. This leads to the 
observed r-process peaks in isotopic abundances Q . 

As already mentioned there are a huge number 
of nuclear reactions that go into modeling a super- 
nova and its aftermath, from the formation of the 
presupernova core, to the high neutron flux and 
rapid capture of neutrons after the collapse. The 
nuclear physics is challenging, in part because of 
the sheer number of reactions as well as the em- 
pirical inaccessibility of some of the unstable inter- 
mediate nuclides. The challenges to theory include 
modeling medium mass nuclei, typically A ~ 50-90 
and at the same time handling of a large number 
of excited states involved. The highly-anticipated 
experiments with unstable nuclei at radioactive ion 
beam accelerator facilities under construction and 
being upgraded around the world |^ will strongly 
test and constrain the models. 

The main task is to find a model with an appro- 
priate number of degrees of freedom. Any model 
with too few microscopic degrees of freedom would 
be inadequate; conversely, a model with too many 
degrees of freedom (such as full-blown shell model 
calculations, which for iron-region nuclides could 
easily require on the order of 10^ basis states) could 
be computationally intractable. 

Fuller, Fowler and Newman [Q (hereafter refered 
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to as FFN) did the most extensive calculations of 
stellar weak rates over a wide range of tempera- 
ture (10^ < T(K) < 10") and density (10 < pY^ 
(gcm~'^) < 10^^). They calculated stellar electron 
and positron emission rates and continuum elec- 
tron and positron capture rates, as well as the as- 
sociated neutrino energy loss rates for 226 nuclei 
with masses between A = 21 and 60. Measured nu- 
clear level information and matrix elements avail- 
able at that time were used and unmeasured ma- 
trix elements for allowed transitions were assigned 
an average value of log ft =5. To complete the 
FFN rate estimate, the Gamow- Teller contribution 
to the rate was parameterized on the basis of the 
independent particle model and supplemented by a 
contribution simulating low-lying transitions. Soon 
it was realized that calculations of presupernova 
evolution generate cores that are so neutron-rich 
{Ye < 0.42) that nuclei more massive than A = 60 
must be considered Q. The FFN rates were then 
updated and extended to heavier nuclei Beta 
decays previously were thought not to play an im- 
portant role in the process of the core collapse and 
were not included in evolutionary calculations (see 
for example 0, 0). It was shown in the work 
of p that beta decays indeed play an important 
role in the evolution of the core below Yg ^ 0.46 
as the core approaches a state of dynamic equilib- 
rium between electron captures and beta decays. 
Motivated by the ideas of these papers || Kar 
et al. calculated f3 decay rates for 11 nuclei 
in the mass range A > 60 using an average beta 
strength function for typical presupernova matter 
density (p = 3 x 10*" to 3 x 10^ g cm~'^) and tem- 
perature (T = (2 to 5) X lO^K). 

First microscopic calculations |l|, |l|, |l|, of 
the terrestrial rates used the pn-QRPA theory, and 
were found to give good agreement with the exper- 
imental data. The pn-QRPA theory is a good com- 
promise for the calculation of weak interaction rates 
in stellar environment between the phenomenologi- 
cal calculations described above and the full-blown 
shell model calculations where the number of basis 
states is too big to handle. The shell model calcu- 
lations can be done on price of truncation of the 
model space but then there is a chance of missing 
certain nuclear effects. An alternate method of so- 
lution, Monte Carlo evaluation of path integrals, 
also known as Monte Carlo shell model ||l^ , allows 
access to larger model spaces but at the cost of 



restricted spectroscopy and a limited class of inter- 
actions. 

This work is the first ever extensive calculation 
of stellar weak rates in the fp- and fpg-sheW nuclei 
ranging from A = 40 to 100 using the pn-QRPA 
theory (see also [|l8|). A total of 619 nuclei were 
considered for the calculation of stellar weak rates. 
These include also proton-rich and neutron-rich nu- 
clei. This work is a continuation of our first pa- 
per for /p-shell nuclei where twelve different weak 
rates had been calculated for each parent nucleus. 
These included e^-capture rates, /3^-decay rates, 
neutrino (anti-neutrino) energy loss rates, energies 
of beta-delayed proton (neutron) and the proba- 
bilities of these beta-delayed particle emission pro- 
cesses. 

Section 2 reviews briefly our formalism of calcu- 
lations which we think is necessary to restate. This 
can help in independent study of this paper. For de- 
tails of formalism we refer to Results, discus- 
sions and comparisons with previous calculations 
are dealt with in Section 3. We finally summarize 
our conclusions in Section 4. 

2 FORMALISM 

The following main assumptions are made: 

(1) Allowed Gamow- Teller and superallowed 
Fermi transitions were taken into account for our 
calculations. Contributions from forbidden transi- 
tions were assumed to be negligible. 

(2) The atoms are totally ionized as a result of 
the high temperature prevailing in the stellar inte- 
rior, and the electrons obey the energy distribution 
function of a Fermi gas. At high temperatures (kT 
> 1 MeV), positrons appear via electron-positron 
pair creation, and the positrons follow the same en- 
ergy distribution function as the electrons. 

(3) Neutrinos and antineutrinos escape freely 
from the interior of the star. Therefore, there are 
no (anti)neutrinos which block the emission of these 
particles in the capture or decay processes. Also, 
(anti)neutrino capture is not taken into account. 

(4) The distortion of electron (positron) wave- 
function due to the Coulomb interaction with a nu- 
cleus is represented by the Fermi function in phase 
space integrals. 

(5) We consider the effect of particle emission 
processes from excited states. All excited states. 
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with energy greater than Sp or Sn, are assumed to 
decay to the ground state by emission of protons 
or neutrons, respectively. Due to uncertainties in 
the calculation of energy levels and the effect of the 
Coulomb barrier for the case of proton emission, it 
is assumed that particles are emitted at excited en- 
ergies 1 MeV higher (called Ecrit throughout this 
section) than the particle decay channel. For the 
parent nuclei, the cut-off excitation energy is set at 
Ecrit- AH excited states in daughter nuclei, with 
excitation energy (Ej < Ecrit) decay directly to 
the ground state through 7 transitions. All excited 
energy states (Ej > Ecrit) decay by emitting par- 
ticles (protons or neutrons) with increasing kinetic 
energies. It is further assumed that either protons 
(if Sp < Sn) or neutrons (if S'„ < Sp) are emitted 
from these excited states. The /? decay of a possible 
isomer is not taken into account. 

2.1 Rate Formulae 

The decay (capture) rate of a transition from the 
ith state of a parent nucleus {Z, N) to the jth state 
of the daughter nucleus {Z ^l, N ±1) is given by Q 



fijiT,p,Ef) 



where {ft)ij is related to the reduced transition 
probability Bij of the nuclear transition by 



^(GT)., = ^^|<j|lE4'^'ll*>l'' (6) 

k 

In Eq. (6), a'' is the spin operator and stands 
for the isospin raising and lowering operator. The 
value of D=6295 s is adopted and the ratio of the 
axial- vector ((7^) to the vector (gy) coupling con- 
stant is taken as 1.254. 

The phase space integral {fij) is an integral over 
total energy. 



fij — 



w^/w^ - l{w„^-w)^F{±Z,w){l-G^)d- 

(7) 

for electron (upper signs) or positron [lower signs) 
emission, or 



fij — 



vj^ — l(w„i + w)'^F{±Z, w)Gzfidw, 

(8) 

for continuum positron {lower signs) or electron 
(upper signs) capture. 

In Eqs. (7) and (8), w is the total kinetic en- 
ergy of the electron including its rest mass, wi is 
the total capture threshold energy (rest-|-kinetic) 
for positron (or electron) capture. One should 
note that if the corresponding electron (or positron) 
emission total energy, Wm, is greater than -1, then 
wi ~ I, and if it is less than or equal to 1, then 
wi =1 Wm |. Wm is the total /3-decay energy. 



(ft),, = D/B,. 



(2) 



md 



Ei — Ej 



(9) 



The D appearing in Eq. (2) is a compound expres- 
sion of physical constants. 



D 



2ln2h\^ 



(3) 



and, 



i?y = B{F),, + {gA/gv)^B{GT),„ (4) 



where B(F) and B(GT) are reduced transition 
probabilities of the Fermi and Gamow- Teller tran- 
sitions respectively. 



B{F),, 



1 



2J, + 1 



\<J\\Y^tl\\^>\\ (5) 

k 

^Throughout subsection 2.1 natural units {fi = c = me = 
1) are adopted, unless otherwise stated, where me is the 
electron mass. 



where rup and Ei are mass and excitation energies 
of the parent nucleus, and rud and Ej of the daugh- 
ter nucleus, respectively. 

6*+ and G- are the positron and electron distri- 
bution functions, respectively. Assuming that the 
electrons are not in a bound state, these are the 
Fermi-Dirac distribution functions. 



G- = [exp{ 



kT 



G4 



[exp{ 



E + 2 + Ef 
kT 



) + iY 



(10) 



(11) 



Here E = {w — 1) is the kinetic energy of the elec- 
trons, Ef is the Fermi energy of the electrons, T is 
the temperature, and k is the Boltzmann constant. 

The number density of electrons associated with 
protons and nuclei is pYgNA, where p is the baryon 
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density, is the ratio of electron number to the 
baryon number, and Na is the Avogadro number. 



pYe 



7 (G_ - G+)p'dp, (12) 
Jo 



where p = {vJ^ — 1)^^^ is the electron or positron 
momentum, and Eq. (12) has the units of moles 
cm^^ . This equation is used for an iterative calcu- 
lation of Fermi energies for selected values of pY^ 
and T. 

In the calculations, the inhibition of the final neu- 
trino phase space is never appreciable enough that 
neutrino (or anti-neutrino) distribution functions 
had to be taken into consideration. F{±Z, w) are 
the Fermi functions and are calculated according to 
the procedure adopted by Gove and Martin | |2C| ]. 

There is a finite probability of occupation of par- 
ent excited states in the stellar environment as a re- 
sult of the high temperature in the interior of mas- 
sive stars. Weak interaction rates then also have a 
finite contribution from these excited states. The 
occupation probability of a state i is calculated on 
the assumption of thermal equilibrium, 



P, 



(2J, + l)exp{-E,/kT) 
E;=ii2J, + l)expi-E,/kT)' 



(13) 



where Ji and Ei are the angular momentum and 
excitation energy of the state i, respectively. 

Unfortunately one cannot calculate the Ji's in 
QRPA theory and hence Eq. (13) is modified as 



P, 



exp{-E,lkT) 
Y.t=i exp{-E,/kT)' 



(14) 



This approximation is a compromise and can be 
justified when one takes into consideration the un- 
certainty in the calculation of Ei which easily over- 
sheds the uncertainty in calculating the values of 
Ji in the above Eq. (13). 

The rate per unit time per nucleus for any weak 
process is then given by 



(15) 



The summation over all initial and final states is 
carried out until satisfactory convergence in the 
rate calculations is achieved. 

The neutrino energy loss rates are calculated us- 
ing the same formalism except that the phase space 



integral is replaced by 

/Win 

(16) 

for electron (upper signs) or positron [lower signs) 
emission, or by 
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l{wm + u')'^i^(±, w)G^dw, 



(17) 

for continuum positron (lower signs) or electron 
(^upper signs) capture. 

The Xij in Eq. (15) in this case is the sum of 
the electron capture and positron decay rates, or 
the sum of the positron capture and electron decay 
rates, for the transition i j. 

The proton energy rate from the daughter nu- 
cleus is calculated, whenever Sp < Sn, by 



A'' — PiXij{Ej ~ Ecrit), 



(18) 



for all Ej > Ecrit, whereas for all Ej < Ecrit one 
calculates the 7 heating rate. 



7 , PiKjEj. 



(19) 



If on the other hand, Sn < Sp, then the neutron 
energy rate from the daughter nucleus is calculated 
using 

X" —\ PiXij{Ej — Ecrit)-, (20) 



for all Ej > Ecrit, and for all Ej < Ecrit the 7 
heating rate is calculated as in Eq. (19). 

The probability of /3-delayed proton (neutron) 
emission is calculated by 



pp{n) _ 



Sij/ Pi^ij' 
Sii Pi^ij 



(21) 



where jf are states in the daughter nucleus for 
which Ej, > Ecrit- In Eqs. [(18)-(21)] A^jf,) is 
the sum of the electron capture and positron de- 
cay rates, or the sum of positron capture and elec- 
tron decay rates, for the transition i — s- j{jf). For 
a detailed formalism of the calculation of excited 
states and the Gamow- Teller transition amplitudes 
between these states as well as the Fermi transition 
amphtudes we refer to [p3. 
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2.2 Incorporation of Experimental 
Data 

The latest experimental excitation energies and log 
ft values are incorporated in our calculations wher- 
ever possible to increase the reliability of our re- 
sults. 

The calculated excitation energies were replaced 
with the measured ones when they were within 0.5 
MeV of each other. The log ft value of this energy 
level was also then replaced by the measured one. 
Very low lying states were inserted in the calcula- 
tions together with their log ft values if the theory 
was missing them. Inverse and mirror transitions 
were also taken into consideration. If there ap- 
peared a level in experimental compilations without 
definite spin and parity assignment, no theoretical 
levels were replaced with experimental levels be- 
yond this excitation energy nor were they inserted 
in the calculations. 

The Q-value of each transition as well as Sp and 
Sn of each nucleus were derived using the experi- 
mental mass compilation of Audi et al. (2l|. For 
nuclei where failed togive mass defects predic- 
tions of MoUer and Nix and Myers and Swiate- 
cki were used to derive the corresponding ener- 
gies necessary for the calculations. The two differ- 
ent theoretical mass formulae used reflect the sen- 
sitivity of stellar rates on input masses. A total of 
31 out of the 209 nuclei (83 out of the total 619 
for all fp- and fpg-sheW nuclei considered in our 
calculations) presented here were calculated using 
theoretical mass formulae. Stellar rates using mass 
formula of Myers and Swiatecki ||2^, for the case 
of neutron-rich and proton-rich nuclei, are given in 
this paper on an abbreviated scale of density and 
marked accordingly (Table B). In |Q the stellar 
rates of neutron-rich nuclei using the mass model 
of MoUer and Nix j22j are presented. The complete 
rate table on magnetic tapes can be requested from 
the authors. 

Table (1) shows the references from which the 
latest experimental excitation energies and log ft 
values are taken (NP Nuclear Physics and NDS 
Nuclear Data Sheets). For those nuclei where 
more than one references were given, the latest data 
was taken into account. 



3 RESULTS AND DISCUS- 
SION 

As stated earlier there are not many calculations 
available for stellar weak rates for nuclei A > 60. 
A compilation of beta-decay rates for 11 nuclei in 
the mass range A > 60 was done by Kar et al. 
]l^ (which we will refer to as KRS in this paper). 
Some assessment of electron capture and (3 decay 
rates for these heavy nuclei was also reported by 
Aufderheide et al. [|[ ^ (which we will refer to as 
ABKSV and AFWH, respectively, in this paper). 
Transitions from only the lowest state of the par- 
ent nucleus having an allowed transition were calcu- 
lated in a very crude way by ABKSV. This was also 
admitted by the authors and comparison of these 
rates with the QRPA calculations clearly reflects 
the situation (see section 3.2). At that time no mi- 
croscopic calculations for /p-shell nuclei were avail- 
able. In their later calculations of electron capture 
and /3 decay, AFWH estimated the /3 decay rates in 
a similar fashion as the electron capture rates, and 
phenomenologically parameterized the position and 
strength of the back resonance. This estimate was 
supplemented by an empirical contribution, placed 
at zero excitation energy, which simulates low-lying 
strength missed by the Gamow- Teller resonance. 
Below a comparison of these earlier compilations 
with this work is presented. 

3.1 Comparison with the calculation 
of KRS. 

KRS constructed a model to calculate the beta- 
decay rates using an average beta strength function 
and an electron-phase space factor evaluated for 
typical presupernova matter density (p = 3 x 10^ 
to 3 X 10^ g cm~'^) and temperature (T = 2 x 10^ 
to 5 X 10^ K). For the Gamow- Teller strength func- 
tion KRS used a sum rule calculated by the spectral 
distribution theory, and the centroid of the distri- 
bution was obtained from the experimental data on 
(p,n) reactions. The width, a, of the Gamow- Teller 
strength function had two parts (ct^ — aj^ + cr^, 
with (Jc = 0.157 Z A~^/'^). The parameter crjv was 
fixed by a best fit to the observed half-lives for the 
free decays of a number of A > 60 nuclei to a global 
average value of 6.3 MeV. 

Fig. 1 shows the comparison of beta-decay rates 
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63CO, Odd-Even 



62MN, Odd-Odd 



BD (sec-^-l) 




Figure 1: Comparison of the QRPA beta-decay (BD) 
rates (this work) with those of [12]. Solid lines represent 
the QRPA beta-decay rates while broken lines represent the 
decay rates of [12]. Log(T) is the log of temperature in units 
of Kelvin and BD represents the rates in units of sec^'^. 
The upper graph is plotted at pYe = 1.6 xlO^ g cm~^ and 
the lower graph is plotted at pYe = 1-6 xlO^ g cm~^. 



65CO, Odd-Even 



BD (sec-"-!) 




Figure 2: Same as Fig. 1 but for P decay of ^^Co. 



BD (sec-^-l) 




Figure 3: Similar to Fig. 1. The upper graph is plotted 
at pYe = 1-3 X lO'' g cm~^ and the lower graph is plotted at 
pYe = 1-3 xIQS g cm-3 for /3 decay of s^Co. 



for ^^IMn. One sees that at tlic lower density and 
log(T) = 2, the QRPA rate (calculated in this 
work) is enhanced by more than an order of mag- 
nitude. KRS considered just the transitions from 
the ground state of ^^IVIn whereas the QRPA rates 
considered Ei up to 5.802 MeV. This is the reason 
that the rates of KRS arc so much suppressed. At 
high temperatures there is a finite probability of oc- 
cupation of parent excited states and KRS did not 
consider the transitions from these excited states. 
At higher density {pYf. = 1.6 xlO^ g cm~'^) the 
respective rates are further suppressed. For nuclei 
with low or intermediate Q- values, the beta-decay 
rates tend to go down at high densities due to the 
blocking of the available phase space of increasingly 
degenerate electrons. At high densities, the contri- 
bution to the rates by the excited states is very 
important, especially for the higher temperatures, 
because such states, due to their higher Q-values, 
are still unblocked. In contrast, at lower densities, 
the differences between the ground-state transition 
rates and the rates inclusive of the excited states 
are not equally large. Figures 2-4 show similar 
comparisons but one notes that the rates of KRS 
get in better and better comparison with the QRPA 
rates. In Fig. 5, eg., for the case of beta-decay of 
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62FE, Even-Even 



69CU, Odd-Even 



BD (sec^-1) 




Figure 4: Same as Fig. 3 but for P decay of ^SCu. 



67NI, Even-Odd 

BD (sec^-1) 




BD (sec'--l) 
0.0004 



0.0003 
0.0002 
0.0001 




Figure 5: Same as Fig. 3 but for (3 decay of ^''Ni. 



BD (sec'^-l) 




Figure 6: Same as Fig. 1 but for /3 decay of ®^Fe. 



^^Ni, the QRPA rate is just about a factor 2 greater 

than that of KRS decay rate (at log(T) = 2.0 and 
pY^ = 1.3 xlO^ g cm-3). At log(T) = 2.0 and 
pYe = 1.3 xlO^ g cm~^ the comparison is again 
fairly good. In case of ^"^Ni, KRS went up to Ei 
of 3.680 MeV as compared to 6.785 MeV consid- 
ered by the QRPA rates. As previously described, 
the parameter (Tjv, m KRS, was fixed by a best fit 
to the observed half-lives for the free decays of a 
number of A > 60 nuclei to a global average value 
of 6.3 McV. For certain cases, KRS rates arc also 
enhanced in comparison to the QRPA rates. This 
seems to be the case for ^^Fe and ^^Cu. Fig. 6 and 
Fig. 7 depict these cases, respectively. In both fig- 
ures one notes that at lower density the beta-decay 
rates of KRS exceed the QRPA rates. At log(T) 
= 2.0 the QRPA rate is suppressed by a factor of 
1.7 for the case of ®^Fe and by a factor of 1.9 for 
^^Cu. However at higher temperatures and densi- 
ties, again, the QRPA rates surpass the beta-decay 
rates of KRS as the authors take a few parent ex- 
cited states into consideration. Fig. 8, Fig. 9 and 
Fig. 10 arc the cases where the rates of KRS arc en- 
hanced at all points of temperature and density. In 
Fig. 10 the abscissa represents the log of beta-decay 
rates so in this extreme case the beta-decay rates 
of KRS are enhanced by a factor of 3.9 at log(T) 
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Figure 7: Same as Fig. 3 but for ^ decay of ^^Cu. Figure 9: Same as Fig. 1 but for /3 decay of ^^Fe. 
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Figure 8: Same as Fig. 1 but for decay of ^^Co. Figure 10: Same as Fig. 3 but for P decay of ^^Cu. 
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64CO, Odd-Odd 
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Figure 11: Same as Fig. 1 but for /3 decay of ''"'Co. 



= 2.0 and pYg = 1.3 x 10^ g cm ^ and by a factor 
of 427 at log(T) = 2.0 and pY^ =1.3 xlO^ gcm-^. 
Fig. 11 represents the special case of ^''Co. In this 
case the strength function of KRS over-estimated 
the terrestrial half-life of ^^Co by an order of mag- 
nitude jl2). They were then forced to change the 
shape of the average beta-strength function. For 
this case the authors changed the value of from 
the global value of 6.3 MeV to 7.5 MeV and also dis- 
torted the Gaussian of the Gamow- Teller strength 
distribution introducing a negative skewness in the 
distribution to bring more strength in the ground- 
state domain. One sees that at lower density their 
Gamow- Teller strength function is still not giving 
strong transitions and the QRPA rate is enhanced 
by a factor of 1.9 at log(T) 2.0 and pY^ = 1.6 
xlO'' g cm^-^. Only at pYe ~ 1.6 xlO^ g cm~^ their 
rates start overtaking the QRPA rates at log(T) > 
3.0. 



3.2 Comparison with the calcula- 
tions of ABKSV and AFWH. 

The electron capture rates of certain isotopes of 
copper and cobalt were calculated considering tran- 
sitions from only the lowest state of the parent 
nucleus having an allowed transition S. The au- 




log(T) 



log(T) 



5> 9.6 9.65 9.7 



Fi gure 12: Comparison of the QRPA electron capture 
(EC) rates (this work) with those of [8]. Solid lines rep- 
resent the QRPA electron capture rates while broken lines 
represent the electron capture rates of [8]. Log(T) is the log 
of temperature in units of Kelvin and log(EC) represents the 
log of electron capture rates in units of sec~^. The upper 



graph is plotted at pY^ = 10^ g 
is plotted at pYe = 10* g cm~'^. 



and the lower graph 



thors used many approximations in calculating the 
Gamow- Teller strength function. A comparison of 
all the electron capture rates of nuclei of A > 60 
(from ABKSV) is presented below. 

Figures 12-14 represent the comparison of elec- 
tron capture rates of ^^Co, ^^Cu and ^''Cu, respec- 
tively. ABKSV underestimated the total transition 
strength which they also pointed out in their paper 
Parent excited states till only 0.5061 MeV were 
considered by ABKSV for the case of ^^Co (in com- 
parison this work considers parent excited states up 
to 7.604 MeV). For ^^^Cu and '^^Cu ABKSV consid- 
ered only transitions from the ground-state. As can 
be easily seen their rates are suppressed since the 
Gamow- Teller strengths from higher parent excited 
states are neglected. For ^^Cu and "^"Cu, ABKSV 
calculated the transition strengths from a method 
given in their Appendix A. They still considered 
transitions from the ground-state but calculated 
quite a high value of the Gamow- Teller strength 
(see Table 2 of ||]). Fig. 15 represents the situa- 
tion for ^^Cu. Electron capture rates for ™Cu were 
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68CU, Odd-Odd 



66CU, Odd-Odd 




9.55 9.6 9.65 9.7 



log(T) 




9.55 9. 



9.65 9.7 



log(T) 



-2.2 
-2 . 4 
-2.6" 



Figure 13: 
66Cu. 



Same as Fig. 12 but for electron capture by 



67CO, Odd-Even 




log(T) 



log(T) 



Figure 14: 

67CU. 



Same as Fig. 12 but for electron capture by 




log(T) 



log(T) 



Figure 15: 

68 Cu. 



9.55 9.6 9.65 9.7 



Same as Fig. 12 but for electron capture by 



not stated as a function of temperature and density 
in their paper and hence no comparison could be 
made for this isotope of copper. 

A search for important weak interaction nuclei in 
presupernova evolution was made by AFWH. Ta- 
bles 4-11 (in H]) rank the most important elec- 
tron capture and beta-decay nuclei where rate cal- 
culations were done for specific values of tempera- 
ture, density and Y^. All of these 160 rates were 
re-calculated using the pn-QRPA theory and the 
comparison of the two calculations are presented in 
Tables (2)-(9). It should be noted that AFWH cal- 
culated the electron capture and /3 decay rates for 
/p-shell nuclei in the mass range 48 to 68. 



4 SUMMARY 

This paper, which is the third in a series starting 
from ||l^, calculates microscopically for the first 
time the weak interaction rates in stellar environ- 
ment for a total of 209 /p-shell nuclei in the mass 
range A = 61 to 80. This region of nuclei is a re- 
gion of notorious complexity and the large number 
of excited states involved makes a full-blown shell 
model calculation a formidable task computation- 
ally. On the other hand various phenomenological 
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parameterizations of the Gamow- Teller strength of 
low-lying states were done previously for an insuf- 
ficient number of heavier nuclei considering only a 
few parent excited states. The pn-QRPA theory 
is a good approach for tackling the problem of mi- 
croscopic calculation of weak interaction rates for 
these nuclei. 

The pn-QRPA theory earlier was quite successful 
in calculation of the terrestrial half-lives of beta- 
minus decay and beta-plus (electron capture ) de- 
cay and good comparison was achieved with the 
experimental data. The model works better with 
increasing neutron excess (see Table B and Fig. 6 of 
). This was a promising feature with respect to 
the prediction of unknown half-lives, implying that 
the predictions were made on the basis of a realistic 
physical model. Our calculations also include the 
weak interaction rates of many neutron-rich nuclei 
which play an important role in the core collapse 
of massive stars. 

We plan in future to employ these microscopic 
weak interaction rates in simulations of stellar and 
galactic evolution processes and nucleosynthesis 
calculations. 
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Table (1): References of experimental data incorporated in this calculation (NP — > Nuclear Physics, 
NDS -» Nuclear Data Sheets). 



Mass Number A 


"r> ^ ^ ^.^ 

Kcicrcnce 


61 


NDb ol,2l\ (1992), NDb 38,463 (1983) 


62 


NDS 60,337 (1990) 


63 


NDS 64,815 (1991), NDb 28,559 (1979) 


64 


NDb 78,395 (1996) 


65 


NDS 69,209 (1993), NDS 47,135 (1986) 


66 


NDS 83,789 (1998) 


01 


inUd D4,o/i3 (^iyyij, iNUb oy,(4i (^iyooj 


68 


NDS 76,343 (1995), NDS 55,1 (1988) 


69 


NDS 58,1 (1989) 


70 


NDS 68,117 (1993), NDS 51,95 (1987) 


71 


NDS 68,579 (1993), NDS 53,1 (1988) 


72 


NDS 73,215 (1994), NDS 56,1 (1989) 


73 


NDS 69,857 (1993), NDS 51,161 (1987) 


74 


NDS 74,529 (1995), NDS 51,225 (1987) 


75 


NDS 60,735 (1990) 


76 


NDS 74,63 (1995), NDS 42,233 (1984) 


77 


NDS 81,417 (1997) 


78 


NDS 63,1 (1991), NDS 33,189 (1981) 


79 


NDS 70,437 (1993), NDS 37,393 (1982) 


80 


NDS 66,623 (1992), NDS 36,127 (1982) 
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Comparison of Calculated Electron Capture Rates with Those of Ref. [9] 

(Table (2)-Table (5)) 

Table (2) 

(Electron Capture Rates: p = 5ME + 07, T9 = 3.40, Ye = 0.47) 



Nuclei 


Rates (from [9]) 


Rates (this work) 


"Co 


3.50E-03 


4.56E-04 


55pg 


1.61E-03 


1.47E-03 


Co 


1.41E-01 


3.99E-02 


54Fe 


3.11E-04 


3.07E-06 


56C0 


7.40E-02 


1.14E-02 


53Mn 


2.48E-03 


9.62E-04 


58Ni 


6.36E-04 


7.31E-05 


59Ni 


4.37E-03 


1.32E-03 


6iCu 


3.93E-01 


5.33E-02 


"Ni 


1.94E-02 


4.76E-02 


58C0 


1.04E-02 


2.23E-04 




5.13E-03 


4.34E-05 


5icr 


2.81E-03 


7.19E-03 


62CU 


1.59E+00 


1.68E-02 


52Mn 


2.85E-02 


2.52E-02 


53pg 


2.04E-02 


5.08E-02 


56Ni 


1.60E-02 


4.83E-03 


s^Cu 


l.OlE+00 


1.42E-01 


60CU 


8.39E-01 


1.29E-01 


60Ni 


1.49E-06 


1.70E-07 



Table (3) 

(Electron Capture Rates: p = 1A5E + 08, T9 = 3.80, = 0.45) 



Nuclei 


Rates (from [9]) 


Rates (this work) 


«"Co 


1.27E-02 


1.74E-05 


59C0 


6.57E-04 


2.09E-04 


s^Mn 


1.57E-02 


3.14E-04 


64Cu 


3.69E-01 


2.35E-02 


61 Ni 


1.20E-03 


3.54E-04 


63CU 


1.86E-02 


1.29E-02 


58C0 


3.07E-02 


1.55E-03 


53Mn 


8.97E-03 


5.64E-03 


55Fg 


6.00E-03 


6.42E-03 


55Mn 


2.25E-05 


9.62E-06 


"Co 


1.29E-02 


2.40E-03 


"Fe 


1.84E-05 


1.05E-05 


62Cu 


4.60E+00 


8.76E-02 


Sly 


2.96E-05 


1.27E-05 


56Mn 


2.56E-04 


2.42E-04 


50y 


2.45E-02 


3.80E-03 


56pg 


1.31E-06 


7.14E-07 


60Ni 


2.74E-05 


3.23E-06 


51 Cr 


9.33E-03 


2.64E-02 


53 cr 


2.46E-06 


8.85E-06 
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Table (4) 

(Electron Capture Rates: p = IME + 09, T9 = 4.93, Yg = 0.43) 



Nuclei 


Rates (trom [9J) 


Rates (this work) 


""Co 


2.31E+00 


1.36E-02 


"''Cu 


6.56E-01 


l.lOE+00 


''"'Co 


3.22E-02 


9.34E-03 


'''^Cu 


2.62E-02 


1.56E-02 




2.90E-02 


2.77E-02 




1.81E-02 


5.86E-04 


'^^Co 


3.35E-01 


3.99E-02 


v^o 




o OOP (T? 


67CU 


2.38E-03 


1.25E-02 


<^5Cu 


5.91E-02 


8.25E-02 


Sly 


9.71E-03 


5.78E-03 


48SC 


8.25E-02 


4.68E-03 


64Cu 


l.lOE+01 


2.51E+00 


63Ni 


3.87E-03 


2.47E-03 


55Mn 


8.73E-03 


5.89E-03 


58Mn 


2.94E-03 


3.23E-04 




4.03E-04 


2.94E-04 


59Fe 


1.83E-04 


5.70E-04 


49xi 


1.70E-02 


7.17E-02 


6lNi 


5.07E-01 


3.46E-02 



Table (5) 

(Eloctroii Capture Rates: /; = -i.OlE + 10. T9 = 7.33, = 0.41) 



Nuclei 


Rates (from [9]) 


Rates (this work) 




1.05E+03 


7.45E+02 


eipe 


1.63E+02 


2.41E+02 


49SC 


5.23E+01 


4.49E+02 


63C0 


1.62E+02 


3.69E+02 




8.36E+02 


3.41E+02 




1.44E+02 


4.23E+02 


64C0 


2.40E+02 


6.29E+02 


55v 


9.23E+01 


5.52E+01 


53y 


1.73E+02 


4.27E+02 


59Mn 


1.41E+02 


2.84E+02 


60Mn 


2.55E+02 


7.19E+02 




6.09E+01 


8.49E+01 


59Fe 


7.20E+02 


2.70E+02 


62C0 


1.44E+03 


5.23E+02 


50Sc 


2.55E+01 


7.39E+02 




6.73E+01 


3.02E+01 


68CU 


2.22E+02 


2.12E+03 


54y 


9.66E+01 


5.11E+02 


56 J, 


3.33E+01 


2.9.5E+01 




1.15E+03 


3.40E+02 
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Comparison of Calculated Beta Decay Rates with Those of Ref. [9] 

(Table (6)-Table (9)) 

Table (6) 

(Beta Decay Rates: p = 5.86E' + 07, TQ = 3.40, = 0.47) 



Nuclei 


Rates (from [9]) 


Rates (this work) 


"Fe 


l.lOE-05 


8.14E-08 


s^Mn 


8.81E-06 


7.17E-10 


58Co 


5.67E-06 


1.09E-10 


53 Cr 


5.93E-06 


6.11E-07 


60Co 


4.66E-03 


1.97E-06 




1.09E-02 


1.86E-05 


55Mn 


3.68E-07 


2.01E-06 


56Fg 


1.19E-10 


1.72E-12 




1.34E-09 


9.93E-12 




8.11E-08 


8.65E-08 




6.95E-03 


3.00E-05 


SOy 


1.89E-05 


3.12E-09 


56 Co 


1.67E-08 


O.OOE+00 


58Fg 


1.09E-07 


1.35E-08 


52Y 


1.60E-02 


6.01E-04 


55pg 


1.30E-10 


2.71E-12 


61 Co 


1.36E-03 


6.55E-04 


63Ni 


3.16E-04 


4.56E-07 


54 Cr 


2.90E-07 


3.36E-08 


s^Cu 


1.34E-04 


2.76E-07 



Table (7) 

(Beta Decay Rates: p = 1A5E + 08, T9 = 3.80, Ye = 0.45) 



Nuclei 


Rates (from [9]) 


Rates (this work) 




1.02E-02 


7.53E-05 


61 Co 


1.90E-03 


3.46E-04 


60Fe 


3.33E-03 


1.93E-05 


56Mn 


7.99E-03 


1.25E-05 


60CO 


4.54E-03 


2.48E-06 


52y 


1.23E-02 


3.24E-04 


55 Cr 


1.52E-03 


3.40E-03 


57Mn 


1.36E-03 


7.99E-03 


53v 


5.60E-03 


5.26E-03 


62C0 


7.87E-02 


6.09E-03 


63Ni 


3.02E-04 


3.65E-07 


53 Cr 


1.57E-05 


7.99E-07 


57Fc 


2.65E-05 


1.26E-07 


54 Cr 


1.61E-06 


8.36E-08 


65Ni 


4.06E-03 


5.15E-04 


58Fe 


4.96E-07 


3.13E-08 


58Mn 


1.71E-01 


9.58E-02 


5lTi 


1.07E-03 


1.45E-03 


63C0 


2.05E-02 


2.19E-01 


6iFe 


8.64E-02 


l.lOE-02 



17 



Table (8) 

(Beta Decay Rates: p = IME + 09, T9 = 4.93, Ye = 0.43) 



Nuclei 


Rates (Irom [9J) 


Rates (this work) 




1.25E-01 


3.49E-03 


°'Cr 


3.31E-01 


4.22E-02 




3.50E-01 


6.32E-01 




6.38E-02 


1.56E-03 




2.01E-02 


1.96E-02 


°'Ni 


1.71E-02 


1.34E-03 




6.02E-01 


6.17E-03 


58T\/r„ 
ivin 


±.Dinj-ui 


o 1 op no 


50Sc 


l.lOE-01 


7.04E-04 




1.85E-01 


1.47E+00 


64C0 


2.29E-01 


3.19E-01 


54y 


1.18E-01 


9.85E-03 




2.70E-03 


2.36E-05 


51xi 


2.13E-03 


6.60E-04 


62C0 


7.16E-02 


4.94E-04 


55 Cr 


3.31E-03 


4.90E-04 




2.84E-03 


1.84E-04 


59Fe 


3.50E-03 


1.84E-04 


eope 


3.63E-04 


3.42E-05 


63Fg 


7.18E-01 


5.98E-02 



Table (9) 

(Bda. Decay Rales: p = -iAYLE + 10. TP = 7.33, Ye = 0.41) 



Nuclei 


Rates (from [9]) 


Rates (this work) 


5^Sc 


1.09E-03 


1.17E-06 


62Mn 


7.78E-03 


4.63E-04 


56v 


2.02E-03 


4.70E-06 


53SC 


1.56E-03 


1.50E-04 


51SC 


6.04E-05 


1.21E-06 


57y 


1.84E-03 


6.33E-05 


60Mn 


4.40E-04 


1.74E-06 


51 Ca 


1.47E-03 


1.26E-05 


55X1 


9.70E-04 


7.64E-05 


68C0 


6.98E-04 


1.20E-03 


61 Mn 


1.86E-04 


5.11E-05 


49Ca 


8.80E-06 


3.23E-07 


50Ca 


4.46E-05 


1.27E-06 


55v 


5.75E-05 


1.20E-05 


59 Cr 


3.50E-04 


3.36E-05 


63pg 


1.17E-04 


6.46E-07 


53X1 


1.53E-05 


2.46E-06 


54X1 


4.64E-05 


3.30E-06 


5"SC 


1.61E-05 


1.16E-08 


"•■Co 


3.43E-(J.5 
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EXPLANATION OF TABLE 



TABLE A. fp-Shell Nuclei Weak Rates in Stellar Matter 



The calculated weak interaction rates [Eqs. (15) and (18)-(20)] are all tabulated in logioA. The 
probabilities of /?-delayed proton (neutron) emission [Eq. (21)] are also tabulated in logarithmic scale. 
These probabilities are calculated only to one significant figure and are given up to three places of 
decimal only for designing purposes. All rates listed for a particular direction concern the parent 
nucleus except for the last two columns which concern the daughter nucleus. For each daughter nucleus, 
either the proton energy rate and probability of /3-delayed proton emission is stated (if Sp < S'„) or the 
neutron energy rate and probability of /3-delayed neutron emission is stated (if Sn < Sp). In the table 
-100 means that the rate (or the probability) is smaller than 10~^°°. 



Q 


Mass of parent minus mass of daughter nucleus 




SnP 


Separation energy of neutron for parent nucleus 




SpP 


Separation energy of proton for parent nucleus 




SnD 


Separation energy of neutron for daughter nucleus 




SpD 


Separation energy of neutron for daughter nucleus 




ADen 


log (pYe) (g.cm""^), where p is the density of the baryon, and Ye 
electron number to the baryon number 


is the ratio of the 


T9 


Temperature in units of 10^ K 




EFermi 


Total Fermi energy of electron and positron, including the rest mass (MeV) 


E+Cap 


Positron capture rate (s~^) 
Electron decay rate (s~^) 




E-Dec 




ANuEn 


Anti-Neutrino energy loss rate (MeV.s"'^) 

Gamma ray heating rate (MeV.s""^), tabulated separately for /?+ 
and /3~ direction 




GamEn 


direction 


E-Cap 


Electron capture rate (s""'^) 
Positron decay rate (s~^) 




E+Dec 




NuEn 


Neutrino energy loss rate (MeV.s"-^) 




ProEn 


Energy rate of /3-delayed proton (MeV.s""'^) 




NeuEn 


Energy rate of /3-delayed neutron (MeV.s~^) 




PPEm 


Probability of /3-delayed proton emission 




PNEm 


Probability of /3-delayed neutron emission 
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61 
23 



V 



38 



61 

24 



CR 



Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-5.026 


1.724 


2.506 


1.463 


0.610 


-1.000 


3.0 


3.0 


-1.971 


1.648 


2.421 


1.362 


0.850 


-0.900 


3.0 


30.0 


2.992 


2.161 


4.193 


2.536 


4.000 


-0.040 


7.0 


1.0 


-10.846 


1.720 


2.499 


1.458 


0.565 


-1.000 


7.0 


3.0 


-3.675 


1.643 


2.414 


1.357 


0.820 


-1.000 


7.0 


30.0 


2.989 


2.161 


4.190 


2.535 


3.997 


-0.040 


11.0 


1.0 


-100 


-58.842 


-59.434 


-60.415 


-84.355 


-20.000 


11.0 


3.0 


-42.152 


-20.376 


-20.503 


-21.224 


-28.351 


-7.000 


11.0 


30.0 


-0.858 


-0.542 


0.515 


-0.191 


0.219 


-0.200 








61 
23 


V38^ 


l\ Tl 39 







Aden 


T9 


E-Cap 


E-I-Dec 


NuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-74.359 


-100 


-74.941 


-75.584 


-90.704 


-20.000 


3.0 


3.0 


-25.580 


-100 


-25.673 


-26.488 


-30.370 


-4.000 


3.0 


30.0 
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EXPLANATION OF TABLE 



TABLE B. fp-Shell Nuclei Weak Rates of Neutron-Rich and Proton-Rich Nuclei in 
Stellar Matter 



The weak interaction rates for neutron-rich and proton-rich nuclei are calculated using the mass formula of 
Myers and Swiatecki ||2^. The explanation of the table is similar to that of Table A. The superscript (MS) on 
each parent nucleus indicates that the mass formula ||2^ is used for the calculation of Q-values and separation 
energies. 
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